SINGULARITIES OF AXISYMMETRIC FREE SURFACE FLOWS 

WITH GRAVITY 
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Abstract. We consider a steady axisymmetric solution of the Euler equations 
for a fluid (incompressible and with zero vorticity) with a free surface, acted 
on only by gravity. We analyze stagnation points as well as points on the axis 
of symmetry. At points on the axis of symmetry which are not stagnation 
points, constant velocity motion is the only blow-up profile consistent with 
the invariant scaling of the equation. This suggests the presence of downward 
pointing cusps at those points. 

At stagnation points on the axis of symmetry, the unique blow-up profile 
consistent with the invariant scaling of the equation is Garabedian's pointed 
bubble solution with water above air. Thus at stagnation points on the axis 
of symmetry with no water above the stagnation point, the invariant scaling 
of the equation cannot be the right scaling. A fine analysis of the blow-up 
velocity yields that in the case that the surface is described by an injective 
curve, the velocity scales almost like V X 2 + Y 2 + Z 2 and is asymptotically 
given by the velocity field 

V(\/X 2 + Y 2 ,Z) = c(—jx 2 + Y 2 ,2Z) 

with a nonzero constant c. 

The last result relies on a frequency formula in combination with a concen- 
tration compactness result for the axially symmetric Euler equations by J.-M. 
Delort. While the concentration compactness result alone does not lead to 
strong convergence in general, we prove the convergence to be strong in our 
application. 



1. Introduction 

Consider the steady axisymmetric Euler equations for a fluid (incompressible and 
with zero vorticity) with a free surface acted on only by gravity. Using cylindrical 
coordinates and the Stokes stream function ip (see for example [HI Exercise 4.18 
(ii)]), we obtain the free boundary problem 

div [ — Vip(xi, X2) ) = in the water phase {-0 > 0} (1-1) 

1 2 

— ^\\/^p(x 1, £2)1 = —X2 on the free surface d{ip > 0}; 
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here the original velocity field 

V(X,Y,Z) = ( - — d 2 ^cos??,-— a 2 ^sini?,— ditp 

\ X\ X\ X\ 

where (X,Y,Z) = (xi cos??, x\ sin??, x 2 )- 

Observe that the positive sign of -0 is chosen just for convenience and that re- 
placing ip by —tp our analysis covers the case of negative ijj as well. 

Note also that the equations above describe apart from a model, where the fluid 
is pumped in or sucked out at a fixed boundary, also the case of a traveling wave 
traveling in the direction of the axis of symmetry; here the equations describe the 
steady flow in the moving frame, so that the original velocity field is 

V(X, Y, Z, t) = V(X, Y. Z- cot) + (0, 0, c ), 

where c is the speed of the traveling wave and 

V(X,Y,Z) = ( - — <9 2 ?/>cosi?,- — d 2 tf)smd, — 

\ Xi Xi Xi 

[T"7] and [19], [18] are excellent reviews on two-dimensional water waves. 



The free boundary problem (1.1) has been studied in [5] where regularity away 
from the degenerate sets {x\ = 0} (the axis of symmetry) and {x 2 = 0} (containing 
all stagnation points) has been shown for minimizers of a certain energy. 

In the present paper we will focus on precisely those two sets and analyze the 
profile of the velocity vector field close to points in those sets. 

Due to the degeneracy of the free boundary condition \Vi/j(xi, 3^2) | 2 = xi 2 X2 at 
points x° = (xIjX®) with a^x!^ = 0, we obtain four invariant scalings 

in the case x\ ^ and x\ 7^ 0, 

in the case x\ and x® = 0, 

in the case x\ — and x\ ^ 0, 

in the case x° = x% = 0. 



ip(x a + 


rx) 


r 




ip(x° + 


rx) 
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rx) 






0(X° + 


rx) 



Note that the velocity (in the moving frame) would scale like 1, |x| 2 , 1, \x\ 2 in the 
respective cases. 

In a first main result we determine the profile of the scaled solution as r — > 



(Proposition 3.10): In the case x\ 7^ and x^^tO the only asymptotics possible is 
constant velocity flow parallel to the free surface. In the case x\ ^ and x 2 = 
the only asymptotics possible is the well-known Stokes corner flow (see [1], |15j . 
|16| . |21|). Due to the perturbed equation the situation is actually not unlike the 
two-dimensional problem in the presence of vorticity (see [20], [5], [6], [7] for two- 
dimensional results in the presence of vorticity) . In the case x° = and x° ^ the 
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only asymptotics possible is constant velocity flow in the gravity direction. This 
suggests the possibility of air cusps pointing in the gravity direction (Figure [l| . 



(0,0) 







FIGURE 1. Dynamics suggested by our analysis 
In the = the only asymptotics possible is the Garabedian pointed 

bubble solution with water above air (cf. [10], Figure [2]). This comes at first as a 
surprise as it means that there is no nontrivial asymptotic profile at all with air 
above water and with the invariant scaling. However there remains at this stage 
the possibility that the solution has a higher growth than that suggested by the 
invariant scaling. 

In Theorem |3.12| we first analyze the possible shapes of the surface close to 
stagnation points and close to points on the axis of symmetry. Assuming that the 
surface is given by an injective curve and assuming also a strict Bernstein inequality 
(corresponding to a Ray leigh- Taylor condition) we obtain the following result: 

In the case x\ ^ and x\ = the only asymptotics possible are the well-known 
Stokes corner (an angle of opening 120° in the direction of the axis of symmetry), 
and a horizontal point. 

In the case x\ = and x® < the only asymptotics possible are cusps in the 
direction of the axis of symmetry. 

In the case x\ = x\ = the only asymptotics possible are the Garabedian pointed 
bubble asymptotics (an angle of opening sa 114.799° with water above air), and a 
horizontal point. 

A fine analysis of the velocity profile in the last case (x\ — x\ — and a horizontal 
point) is no mean feat, and we confine ourselves to the case of air above water. Here 
we prove (Theorem 7.1l that the velocity scales almost like \JX 2 +Y 2 + Z 2 and is 
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Figure 2. Garabedian pointed bubble asymptotics 



asymptotically given by the velocity field 

V(VX 2 +Y 2 , Z) = c(-Vx 2 + Y 2 , 2Z), 
where c is a nonzero constant (Figure k3| . 



(0,0) 




FIGURE 3. Dynamics suggested by our analysis 
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The proofs rely on a monotonicity formula as well as a frequency formula for 
the axisymmetric problem; as remarked in |21| . it is for certain semilinear problems 
possible to derive on the set of highest density not a perturbation of Almgren's 
frequency formula (see [T], [T3], [H], but a true nonlinear frequency formula. 

Here we extend the formula of |21| to the axisymmetric case. In combination with 
a concentration compactness result for the axially symmetric Euler equations by 
J.-M. Delort jS], this leads to the already mentioned profile for the velocity vector 
field. Note that while the concentration compactness result alone does not lead to 
strong convergence in general, we prove the convergence to the limiting velocity 
vector field to be strong in our application. 

2. Notation 

We will use coordinates (X, Y, Z) in the physical space R 3 together with partial 
derivatives dx,dy,dz as well as two-dimensional coordinates x — {x\ 1 X2) together 
with partial derivatives 81,82- Sometimes we are going to to use cylindrical coor- 
dinates (X, Y, Z) — (xi cost?, Xi sin??, x%). We denote by x ■ y the Euclidean inner 
product in R™ x R™, by |a;| the Euclidean norm in R", by B r (x°) := {x e R™ : 
\x — x°\ < r} the ball of center x° and radius r, by B+(x°) := {x € R" : x\ > 
and I a: - x°\ < r}, by 8B+(x°) := {x 6 R" : 11 > and \x - x°\ = r} and 
R" := {(xx, . . . , x n ) : xi > 0} the positive parts. Note that 8B^(x°) is not the 
topological boundary of B^(x°) and that B+(x°) is not necessarily a half ball. 

We will use the notation B r for B r (0) as well as B+ for i3+(0), and denote by 
UI2 the 2-dimensional volume of B\. 

We will use the weighted LP space 

L^(R+) := {v measurable : [ — \v\ p dx < +00} 

1 

with norm ||/||lp,(r^) = ^/r,2 -^\v\ p dx^j r , the weighted Sobolev space 

W^' P (R^) :={v e L p w (R+) : all weak partial derivatives of v 
are contained in L p (R+)} 
as well as the local spaces 

L v w ;oc (R^_) := {v measurable : / — \v\ p dx < +00 for each R g (0, +00)} 

Jb+ x i 

and 

W w ' p 1oc (R?l_) :={v measurable : v e L p w (oc (R^_) and all weak partial derivatives 
of v are contained in (B^) for each R g (0, +00)}. 

We denote by xa the characteristic function of a set A. For any real number a, 
the notation a + stands for max(a, 0) and a~ stands for min(a,0). Also, C n shall 
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denote the n-dimensional Lebesgue measure and H s the s-dimensional Hausdorff 
measure. By v we will always refer to the outer normal on a given surface. We will 
use functions of bounded variation BV(U). i.e. functions / e ^(U) for which the 
distributional derivative is a vector-valued Radon measure. Here |V/| denotes the 
total variation measure. Note that for a smooth open set E C R 2 , |Vx_e| coincides 
with the surface measure on dE. We will also use the reduced boundary d re dE. 



3. Notion of solution and monotonicity formula 

Let n be a bounded domain contained in {(xi,x 2 ) '■ x\ > 0}, in which to 
consider the combined problem for fluid and air. We study solutions u, in a sense 
to be specified, of the problem 



div ( — Vu) 

Xl 



8 





du \ 


a 




du \ 


(xi 


dxi) 


dx 2 




dx 2 J 



in fin {w > 0}, (3.1) 



1 9 

-2|Vu| = x 2 onnnd{u> 0}. 

X-i 



Note that, compared to the Introduction, we have switched notation from tp to u, 
and we have "reflected" the problem at the hyperplane {x 2 = 0}. Since our results 
are completely local, we do not specify boundary conditions on dQ. 



We begin by introducing our notion of a variational solution of problem (3.1 ) 



Definition 3.1 (Variational Solution). We define u € W^ oc (Q) to be a variational 



solution of Q if u G C°(fi) H C 2 {tt n {u > 0}), Vu/xi eC 1 (Cln{u> 0}), it = 
on {xi = 0} (motivated by the fact that the velocity on the axis orthogonal to the 
axis direction should be zero), u > in 17, and the first variation with respect to 
domain variations of the functional 

■ r \-<'2\\r n] 



J(v) 



dx 



vanishes at v = u, i.e. 



0= - — J{u{x + ecj)(x))% =0 



1 



XiX 2 X{u>Q} 



div 



2 — VuD<j)Vu 

Xi 



for any cf> 



■ I l Vtt l + 2; 2X{n>0}J 01 + ^lX{u>0}^2 

2 ) € C^(Q; R 2 ) such that X = on {x x = 0}. 



dx 



A proof of the just mentioned first variation formula can be found in |14| Section 
3.2] . An integration by parts shows that u satisfies on smooth parts of free boundary 
d{u > 0} in {xix 2 ^ 0} the free boundary condition 

1 



xi 



Vtt = xix 2 . 
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Theorem 3.2 (Monotonicity Formula). Let u be a variational solution of (3.1), 
let x° € Q, and let 5 := dist(a; , dCl)/2. Let, for any r £ {0,5), 

I{r)= ( ( — \Vu\ 2 + x 1 x 2X{u >q} ) dx, (3.2) 

JB+(x°) \ x 1 J 



B + (xO) 
f 

dB+(x°) x l 



Jin = I ^-u^dU 1 , (3.3) 

M mt (r) = r- 2 I(r) - r~ 3 J(r), (3.4) 

M X2 ( r ) = r -3/( r ) _ 3 r - 4 J(r), (3.5) 

M Xl (r) = r- 3 I{r) - 2r~ 4 J(r), (3.6) 

M XlX2 (r) = r- 4 I(r) - ^r~ 5 J(r). (3.7) 
Then, for a.e. r € (0,5), 

(M int (r))' =2r~ 2 [ — (vu ■ v - -)* dH 1 (3.8) 

JdB + (x«) x l V rJ 

+ r- 3 [ (- Xl ~ Xl \Vu\ 2 + [(a?i - x?)a; 2 + (ar 2 - a:§)ai] X{«>o}) dx 

JB+(x°) \ x 1 J 



- r 



lB+(x°) 

I Xl ~ x Wdu\ 



OB 



+ {x°) ( x lY 



In the case x\ = 0, 



1 / 3u" 2 



(M X2 (r))' = 2r~ 3 / — (Vu-v---) dH 1 (3.9) 

JdB+(x°) x x\ 2rJ 

+ r~ 4 / (- — Xl \Vu\ 2 + (xi - x° 1 )x 2 X{u>o}] dx 
Jb+(x°) V x i J 

3 5 r xi — Xi 2 



J+(x°) ( x l) 



In the case X® = 0, 

(M Xl (r))' = 2r- 3 / — f V« • v - 2-) ' dH 1 (3.10) 

JdB+(x°) x l V r/ 

+ r~ 4 (x 2 - x°)xiX{ u >o} dx. 

JB+(x°) 

Last, in the case x\ = x\ = 0, 
(M XlX2 (r))' = 2r- 4 f — (vu-v- 5u \ dH 1 . (3.11) 



Remark 3.3. (i) The integrand in the first integral on the right-hand side of (3.8) 
is a scalar multiple of (Vh(i) • (x — x a ) — u(x)) 2 , and therefore vanishes if and only 
if u is a homogeneous function of degree 1 with respect to a; . 
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(ii) The integrand in the first integral on the right-hand side of (3.9 1 is a scalar 



multiple of (Vti(s) • (x — x°) — |u(x)) 2 , and therefore vanishes if and only if u is a 
homogeneous function of degree 3/2 with respect to x°. 



(iii) The integrand in the first integral on the right-hand side of (3.9 1 is a scalar 
multiple of (Vit(i) • (x — x°) — 2u(x)) 2 , and therefore vanishes if and only if u is a 
homogeneous function of degree 2 with respect to x° . 



(iv) The integrand in the first integral on the right-hand side of ( 3.9 ) is a scalar mul 



tiple of (yu(x)-x— ^u(x)) 2 , and therefore vanishes if and only if u is a homogeneous 
function of degree 5/2. 

Proof. First, for each u € W w ' ; oc (R 2 ), each a G R and a.e. r £ (0,S) we obtain, 
setting w r (x) := u(x° + rx), 



r ! — u 2 dH 1 ) = — I /■' 

dr V JdB+(x°) x l J dr \ 

1 



as+ x i + ra; i 



w 2 ^ 1 (3.12) 



x x 



{a + n-iy - 1 — u'dH 1 -r a+n - 1 . 

JdB+{x°) x i JdB+ i x i + rxi) 

2 

-w r Vu(x + rx) ■ x dH 1 



w 2 dH 1 



2 r 



a+n—1 



+ xl + rx\ 



dBT x l 



(a + n — 1) 



,a-l 



1 



u 2 dU X -r a - 1 



dB + (x°) x l 



3+(x°) ( x l) 



X l X l „ 2 7,1/1 

u dH 



I — uVh • v dH ■ 

3+(x°) x l 



Suppose now that u is a variational solution of (|3.1|). For small positive k and 

x — 



i] K (t) : = max(0, min(l, L -^)), we take after approximation 4> K (x) := r) K (\x — x |)( 
x ) as a test function in the definition of a variational solution, obtaining 



„0\2 1 



(dm) 2 ^Qx - x°\) + r,' K (\x - ^°l) ( ^:^| 2 ) 
+ (d 2 u) 2 (v K (\x - x°\) + r,' K (\x - A) { *lZf* 

+ (d 1 u)(d 2 u)2r 1 '(\x - x°\) {Xl - ^)(^ -^) 



+ 



/ f " ^-2^|V«| VI* - Z°|) + [(.X! - X?)^ 

+ (x 2 - x%)xi] X{u>o}Vn(\x - x°\)j dx. 



dx 
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Passing to the limit as n — > 0, we obtain for a.e. r £ (0,(5), 



0=2/ xxX2X{u>o}dx-r [ ( — |Vw| 2 + x 1 x 2 X{u>o} ) dH 1 (3.13) 

JB+(x°) JdB+(x°) \ x l J 

2r I — {Vu-vf dV. 1 

JdB+(x°) x l 



X\ — X 



.0 

-|Vu| 2 + [(xi - x?)x2 + {x 2 - xl)xi]x{ u >o} 



x\ 

Observe that letting e — > in 
/ — V«- Vmax(u-e, 0) 1+e dx = 

JB+(x°) x 1 

for a.e. r G (0,<5), we obtain the integration by parts formula 



dx. 



f — max(u - e,0) 1+e Vu • dU 1 

dB+(x°) x l 



for a.e. r G (0, <5). 
Note that 



/ — \Vu\ 2 dx= [ —uVu-vdV 1 

JB+{x°) x 1 JdB+(x°) x l 



(3.14) 



(r- 2 /(r))'=- 



= -2r 3 / ( — |Vu| 2 + xix 2 X{„>o} ) 

_-2 



1 

ldB+{x°) \ x l 



|Vu| 2 +x 1 x 2 X{u>o} ) rf'H 1 



so that by pTl3j ) and pil l 

(r- 2 /(r))' = r~ 3 



[2r / — {Wu-vf dU 1 -2 / —uWu-vdU 1 

\ JdB+(x°) x l JdB+(x«) x l 

( Xi - X° \ 

2 - ± W U \ 2 + [i X l ~ X l) X 2 + ( X 2 ~ x 2) x l] X{«>0} 

.0) \ X l J 



(3.15) 

dx , 



Jb+(x°) 

Combining (3.151 and (3.121 with a = —3 yields (3.8 1 
Moreover, 

' 1 ( — |Vm| 2 +XiX 2 X{u>0} ] 



(r~ 3 I(r))' = -3r~ 



— \Vll\ 2 + XiX 2 X{u>0} 
dB+(x°) \ x l 



dx (3.16) 
dH 1 . 



In the case x 2 = we obtain from (3.161, using (3.131 and (3.14), that 



(r- 3 I(r)Y = r 



l (2r / — (Vu • i>) 2 dft 1 - 3 / 

\ JdB+(x°) x l JdB+(x 



— uVh ■ v dH 1 

■o) Xi 

(3.17) 



/B+(x°) 



x\ 



|Vu| 2 + (xi - x i)x 2 X{u>o}^ 



dx 



10 



E. VARVARUCA AND G.S. WEISS 



Combining (3.171 and (3.12) with a = —4 yields (3.9). On the other hand, in the 
case x\ = we obtain from (3.16), using (3.13) and (3.141, that 



(r- 3 I{r)Y = r- 4 2r 



\2r [ — {\7u- v) 2 dH 1 - 4 [ —u^/u-vdl-L 1 

\ JdB+(x°) x l JdB+(x°) x l 



(3.18) 



/ (%2 ~ A) x iX{u>o} dx , 

JB+(x°) I 



Combining (3.181 and (3.121 with a = -4 yields (3.101. 



Last, in the case X® = x\ = 0, since 



(r- 4 /(r))' = -4r- 



B + (0) \ x l 



1 



Vu| + xix 2 X{u>o} dx 



as+(o) 



^|Vu| 2 +x 1 x 2 X{«>o}^ dU 1 , 



we obtain from (3.13) and (3.141 that 



(r- 4 /(r))' = r~ 



\2r — (Vu-t^tfH 1 -5 / —u{Vu-v)dU 1 

\ JdB+(0) x l JdB + (0) x l , 



Combining ( |3.19| and \i.\2\ with a = -5 yields ( |3.11 l. 

Lemma 3.4 (Bernstein estimate). In {u > 0}, £/ie solution satisfies 



(3.19) 
□ 



A 



Proof. Direct calculation. 



|Vm| 

Xl 



XiX 2 



= 2E 



«,j=l 



□ 



Remark 3.5. Constructing barrier solutions it is therefore possible to verify — 
x\X2 < for certain domains C \xi > 0}, certain Dirichlet boundary data and the 
minimal solution u (cf. [22 ). 



Definition 3.6 (Weak Solution). We define u € W^ 2 oc {VL) to be a weafc solution 
of (3.1l if the following are satisfied: u is a variational solution of (3.1) and the 
topological free boundary d{u > 0} n fl° n {x 2 7^ 0} is locally a C 2,Q -surface. 

Remark 3.7. (i) It follows that in Sl° n {a;2 ^ 0} the solution is a classical solution 
of ([31]). It follows also that d{u > 0} C {x 2 > 0}. 
(ii) For any weak solution u of (3.1 1 such that 

< Cxila^l locally in fi, 
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it is a variational solution of (3.1 1, X{u>o} is locally in Q° n {x2 > 0} a function of 
bounded variation, and the total variation measure |Vx{ u >o}l satisfies 



/ v4 d\Vx{ u >o}\ < Ci 

JB+(x°) 



r 2 , x\ = 



i v z i n.i»^wji — —j- i / (J ^ n 

)B+{x°) v [ rV^'^ > 

for all B+(x°) CC ft. The reason is that, integrating by parts, 

0=/ div(^V U )</ ^d^ 



B+(z )n{u>0} x l J dB+ (x°)n{«>0} 



B r (i°)n{xi=o} 35 1 Js+(cc )na red {M>o} x i 



\y0B+(x°) JB r (i»)n{n=0} / 

Jx+dU 1 . 

; + (x°)n3 r „ d {u>0} 

Lemma 3.8. Let u be a variational solution of (3.1 \ and suppose that 

|Vu| 2 

< C^CiIxqI locally in f2. 

Xl 

Then: 

(i) The limit M int (Q+) = ]im r ->o+ M int (r) exists and is finite. If xl = 0, then 
the limit M X2 (0+) = lim,._>o+ M X2 (r) exists and is finite. If x\ — 0, then the limit 
M Xl (0+) = lim J ,_ i .o+ M Xl (r) exists and is finite. If x® — x^ — 0, then the limit 
M XlX2 (0+) = lim r _i.o+ M XlX2 (r) exists and is finite. 

(Hi) Let xl > 0, x® > and < r m — > 0+ as m — > oo be a sequence such that 
the blow-up sequence 

u m (x) := u(x° + r m x)/r m (3.20) 

converges weakly in W, ' c (R 2 ) to a blow-up limit u Q . Then Uq is a homogeneous 
function of degree 1, i.e. uq(Xx) = Xuq(x). 

Let x\ = and let < r m — > 0+ as m — > oo be a sequence such that the blow-up 
sequence 

u m (x) := u(x° + r m x)/r m i (3-21) 

converges weakly in W^ c (R 2 ) to a blow-up limit uq. Then uq is a homogeneous 
function of degree 3/2. 

Let Xi=0 and let < r m — > 0+ as m — > oo be a sequence such that the blow-up 
sequence 

u m (x) := u(x a + r m x)/r m 2 (3.22) 

converges weakly in W w ' loc (R+) to a blow-up limit uq. Then uq is a homogeneous 
function of degree 2. 
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Let x\ = x\ = and let < r m — > 0+ as m — > oo be a sequence such that the 
blow-up sequence 

u m (x) := u(x° + r m x)/r m i (3.23) 

converges weakly in W^'| oc (R^_) to a blow-up limit Uq. Then uq is a homogeneous 
function of degree 5/2. 

(Hi) Let u m be one of the converging sequences in (ii). Then u m converges 
strongly in W v J ]loc (Bj+) (strongly in W x ^ c (R 2 ) in the cases where x\ > 0). 

(iv) If x® > and x° 0, then 



M mt (0+) = x»x» lim r~ 2 / X{u>0} dx. 
r ^°+ JB+(x°) 

Moreover, M mt (0+) = implies that Uq — in R 2 for each blow-up limit uq of 
u m (x) = u(x° + r m x)/r m . 
If x\ > and x\ = 0, then 



M x *(0+)=x° 1 lim r 3 / x 2X{u >o}dx. 

r ^°+ Jb + (xO) 



Ifx\ = Q and x\ ^ 0, then 



M Xl (0+)=x IJ 2 lim r d / x lX{u >o} dx. 



Moreover, M Xl (0+) = implies that uq — in R 2 ^ for each blow-up limit uq of 
u m (x) = u(x° + r m x)/r m 2 . 
If x\ = x° 2 = 0, then 

M x ^(0+)= lim r- A [ Xl x 2 X{u>o}dx. 

Proof, (i) follows from the assumption 

|Vu| 2 

< Cxi|x2| locally in fl 

Xi 

together with Theorem |3.2| 

(ii): For each < a < oo the sequence u m is in each case by assumption 
bounded in C 0,1 (B+) (bounded in C 0,1 (B a ) in the case that x\ > 0). For any 
< r < a < oo, we write the identities (3.8|, (3.9), (3.8), (3.11 1 in integral form as 



r 2 I — (Vu-u- - ) 2 dU x dr 

3+(x°) £i v rJ 



x°) x l V 

= M m *(cr) - M mt {r) - J K mt {r) dr in the case x\ > and x\ > 0, (3.24) 
[ r~ 3 [ — (vu-v~ 3U ] dU X dr 

Jt JdB+(x°) %l V 2 r J 

= M X2 (a) - M X2 (r) - J K X2 (r) dr in the case x\ > and x\ = 0, (3.25) 
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1 

dB+(x°) x l 



Vii-iz-2-) dU l dr 



M Xl (a) - M Xl (t) - J K Xl (r) dr in the case x? = and x% > 0, 



dB+(x°) x l V 2 r J 

M XlX2 (a) - M XiX2 {t) dr in the case x° = x\ = 0; 



(3.26) 



(3.27) 



here K in \K X2 and # Xl are defined by fl3.24| , ( |3.25| and ( |3.26| , and they are all 
integrable. 



It follows by rescaling in ( 3.24 1-( 3.27 1 that 



2 / \x\ 3 — (Vu TO (s) • x — u m (x)) 2 dx 

I bao)\b t (o) Xl 

rint ( „ _n T^A'iUtf^, ^\ i / I Tsint 



< M mt (r m a) - M mt {r m r) + \K lnt (r)\ dr -> as to -> oo, 

in the case x^ > and x" > 0, 

2/ |x|~ 5 — ( Vw m (x) • x - ^u m (x) J dx 

's„(o)\s r (o) x i V 1 J 

<M X2 {r m a)- M X2 (r m T)+ \ K X2 (r) | dr -> asm^oo, 

in the case x^ > and x% = 0, 

2 / M~ 5 — (V u m (x) ■ x — 2u m (x)) 2 dx 

/b+(o)\b+(o) 



< M X1 {r m a) - M X1 (r m r) + \K Xl (r)\dr^0 as m oo, 
in the case x^ 1 = and x\ > 0, 

2/ M~ 6 — ( Vu m (x) • x - -u m {x) ) dx 

Jb+(o)\b+(o) x i V 1 J 

< M XlX2 (r m cr) - M I1I! (r m r) ^ as m -> oo 
in the case x\ = x% = 0, 

which yields the desired homogeneity of Uq. 

(iii): In order to show strong convergence of Vu m , it is in view of the weak 
L^-convergence of Vit m sufficient to prove convergence of the Ll^-norm. 

Let 8 :— dist(x°, dSl)/2. Then, for each m, u m is a variational solution of 



\(x u + r m x)i 



in B§/ rm fl {it m > 0} in the case x\ > 0, 
in Bf^ r n {u m > 0} in the case x\ = 0. 



(3.28) 
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Since u m converges to uq locally uniformly, it follows from (3.28 1 that uq is harmonic 
in {uq > 0} in the case x® > and a solution of the equation 



div Vu j = 

in the case x\ — 0. Also, using the uniform convergence, the continuity of u and 
its solution property in {uo > 0} we obtain as in the proof of (3.141 that 

+ / -\\^u m \ 2 T]dx 
JR 2 x i 

= [ , n , 1 T-\^u m fr\dx = - [ u m 1 —\7u m -\7r)dx 

J R 2{x u + r m x) 1 J R 2 (xu + r m x)i 

— > — [ uq— frVtJo • Vrjdx = — ?t / I Vito| 2 ?7 dx in the case x? > and that 

/ — |Vu m | 2 r7 g?x = — / u m — Vu m -Vr]dx 
Jul X\ J R 2 Xi 

— > — [ uo — Vuo • \7r]dx — [ — |Vuo| 2 f?<ix in the case x? = 
Jul xi ' Jul xi 

as to — > oo. It therefore follows that Vu m converges strongly in (and in L 2 if 
a;? > 0) to Vuo as m — > oo. 



(iv): Let us take a sequence r TO — > 0+ such that u TO defined in ( 3.20 )-( 3.23) 

1,2 / R 2 
tu,locV xv + 



^' C 2 (R 2 ) in the cas° ■•" 
function uq. Using (iii) and the homogeneity of uq, we obtain that 



converges weakly in W^ C (K%) (weakly in W^' c f(R 2 ) in the case x\ > 0) to 



a 



lim M int (r m ) = ~ ( / \Vu \ 2 dx- f uldli 1 
x\ \J Bl J dBl 



lim r 2 

r->0+ 



£i£2X{«>o} 



x?x!^ lim r 2 / X{«>o} dx, 

\( [ |Vu | 2 dx- 
x i \Jb 1 



lim M X2 (r m ) 



uldU 1 



lim r 3 

r-»0+ 



B+(x° 

= x° lim r~ 3 / £2X{n>o} 
f - y0 + Jb+{x°) 



XlX 2 X{u>0} 



dx 



1 



1 



lim M Xl (r m ) — I —\Vu { 2 dx-2 I —u^dH 1 



dB+ x l 



lim r 3 

r->0+ 



a;i^2X{«>0} 



fix 
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= A !im r 3 / xiX{ u >o}dx 



B + (x°) 

lim M I1IJ (r m )= / — \\7u \ 2 dx- 5 [ —u 2 d% 1 
'™^ 00 Jb+ xi 2 7as+ &i 



r->0+ 



lim r / XiX 2 X{«>0} dx 



B + (xO) 



r->0+ 



lim r / xix 2 X{u>Q}dx 



B+(x°) 



In the case a:§ > 0, M" l *(0+) > and M Xl (0+) > 0, and equality implies that u m 
converges to in measure in R/j_. □ 

The next lemma will be useful in the characterization of blow-up limits in Propo- 
sition mini 

Lemma 3.9. The Legendre function y = P3/2 satisfies 

y'(x) 

x 1— > — r is strictly increasing on (—1,1). 

y'{-x) 

Proof. It suffices to prove that 

y"(x)y'(-x) + y"(-x)y'(x) > in (-1, 1). 
Using the differential equation 

(l- x 2) y »( x )-2xy'(x) + ~y(x) = 0, 

we obtain 

y"(x)y'(-x) + y"(-x)y'(x) = ----—^(y(x)y'(-x) + y(-x)y'(x)). 

4 1-r 

Therefore it is sufficient to prove that f(x) = y(x)y'(—x) + y(— x)y'(x) < in 
(—1, 1). As f(x) — > —00 for |x| —> 1, must have a maximum point in (— 1, 1). At 
the maximum point, 

2r 

= f'(x) = ^y(x)y'(-x) + y(-x)y'(x)), 

1 — ar 

implying that x = and that 

max / = 2v(0y(0) = 3 ^-^- P^O) 

r(-l)r(|)r(i)r(|) < 



(see http://functions.wolfram.com/07.07.20.0006.01 , 



http://functions.wolfram.com/07.07.03.0001.01). □ 
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Proposition 3.10 (Characterization of blow-up Limits). Let u be a variational 

at 

< Cxi\x2\ locally in i7, 



solution o/ (3.1 1 , and suppose that 

|Vu' 2 



Xl 

and that 



i J x t d \^X{u>o}\ < Ci 

JB+(x°) v 



T2 , x\ = 



>B+(x°) v [ r^^, a;^ > 

for all sufficiently small r > 0. 
TTien i/ie following hold: 

(i) In the case xi > and cc!] > 0, the only possible blow-up limits of u m (x) — 
u(x° + r m x)/r m are 

u (x) — x\\Jx\vsiayi(x ■ e, 0) and Uq(x) = j\x ■ e|, 

where e is a unit vector and j is a nonnegative constant. In the case Uq(x) — 
x l \/%2 max ( 2; ' e i0), the corresponding density is M mt (0+) = xlx^^/^t in the 
case uq(x) = j\x ■ e| with 7 > the density is M mt (0+) — xIx^uj^, while in the 
case uq = the density has possible values M mt (0+) £ {0, x^x®^}- 
fiij In the case xi > and £2 = 0, the only possible blow-up limits are 

y/2x° 3 

u Q {p sin 9, p cos 9) = -^p 3 ' 2 cos(-9)x{(p S - m e,pco S e):-7r/3<e<^/3}, 
with corresponding density 

M X2 {Q+)=X^ I X 2 X{(p S me,pco S 0):-7r/3<e<iT/3} dx , 



IB, 

and uq(x) = 0, with possible values of the density 

M X2 (0+) e 1x1 { x%dx,x\ I x^dx,0\. 



(Hi) In the case xi = and x^ > 0, the only possible blow-up limits are 

Uo(x) = yx\ 

with 7 a nonnegative constant and corresponding density 

M Xl (Q+) =x% [ x x dx, 

and Uq(x) = 0, with possible values of the density 

M Xl (0+) e I x\ f xidx,0}. 



(iv) In the case xi — x® = 0, the only possible blow-up limits are 

u (p sin 9, p cos 9) = r% Ue(9) 
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with corresponding density 

M X1X2 {0+) = f Xl x 2 dx, 

J B+ n{(p sine, p cos 0):P^, 2 (- cos0)<O} 

where P3/2 is the Legendre function and Ui is a unique function which is positive 
in H {Py 2 (— cos8) < 0} (an angle of m 114.799° in the positive X2~direction) 
and zero else, and Uq(x) = 0, with possible values of the density 

M XlX2 (0+) £ |y ^xix^dx, J ^XiXz dx,0} . 



For Ui we have the relations 

^U e (8) = c sin 2 6P; /2 (cos6) : Ue(0) = Co|sin0P 3/2 (cos0) 
with a unique positive constant Cq. 

Proof. Consider a blow-up sequence u m as in Lemma |3.8| where r m — > 0+, with 
blow-up limit uq. Because of the strong convergence of Vw m to Vito m L 2 and the 
compact embedding from BV into L 1 , uq is a homogeneous solution of 

0=/ -4f|Vu | 2 div (j) -2Vu Q DfaVu \ dx + x\x% [ xodiv fadx (3.29) 

JR.* x\\ J J R 2 

in the case x\ > and x% > 0, 

0=y (jVw | 2 div (j) - 2Vito-D0V«o) rfx + / ^ir^Xodiv </> + o^Xo^) dx 

(3.30) 

in the case x\ > and a;^ = 0, 
°= / — f|Vu | 2 div 4>- — \Vu \ 2 fa - 2Vm Oj D0Vu o N ) dx (3.31) 

in the case x\ — and £2 > 0j 
0=/ — f|Vu | 2 div (j)- — \\7u \ 2 fa - 2\7u Dfa¥u ) (3.32) 

(xix 2 Xodiv 4> 4- x 2 Xo</»i + ZiXo^J dx 

in the case = a;" = 0; 

the formulas are valid for every <p = (fa, fa) £ Cq(R 2 ; R 2 ) in the case x® > and 
for every </> = (fa, fa) £ Cq(R 2 ;R 2 ) such that fa = on {xi = 0} in the case 
Xi = 0. Moreover \o is the strong Lj^-limit of X{u m >o} along a subsequence. The 
values of the function xo are almost everywhere in {0, 1}, and the locally uniform 
convergence of u m to uq implies that xo — 1 m { u o > 0}- Moreover xo is constant 
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in each connected component of {uq = 0}° \ {x 2 = 0}. In the case uo = 0, (3.29)- 
(3.321 show that xo is constant in {x 2 ^ 0} in the cases ( |3.30 I and (3.321 and that 
Xo is constant in the cases ( 3.29| and (3.31 1. Its value may be either or 1. 

Let z be an arbitrary point in 8{uq = 0} \ {0}. Consider first the case when 
Bg(z) n {uq > 0} has exactly one connected component. Note that the normal to 
8{uq = 0} has the constant value v(z) in Bs(z) for some S > 0. Plugging in <j>(x) :— 



r}{x)v{z) into ( 3.29 1-( 3.32 1 , where 77 € Cq(Bs(z)) is arbitrary, and integrating by 
parts, it follows that 



= 



f — ^IV^qI 2 + x\x%{1 - xo^vdH 1 
u >a} \ x i J 

in the case x\ > and x% > 0, 

/ (- — IV^ + ^I-Xo))^ 1 

/9{u >0} V X l J 

in the case x\ > and x® = 0, 

/ (- — \\7u \ 2 + x lX ° 2 (l -Xo^vdH 1 

Jd{u o >0} \ X l J 

in the case x\ — and x\ > 0, 
1 



0=/ \Vu \ 2 + x lX2 (l - xo) ) vdn 1 

Jd{u >0} V x i 

in the case x\ — x® = 0. 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



Here xo denotes the constant value of xo m { u o = 0}°. Note that by Hopf's 
principle, Vu • v 7^ on B$(z) n d{u n > 0}. In all cases it follows therefore that 
Xo 7^ 1, and hence necessarily xo = 0. We deduce from ( 3.33 1-( 3.36 ) that 



|Vu | 2 =(a;i) 2 a;2 on d{u Q > 0} 

in the case x\ > and x\ > 0, 
|Vu | 2 =( x ifx2 on d{u > 0} 

in the case X® > and x® = 0, 
|Vu | 2 =x 2 X2 on d{u > 0} 

in the case x\ = and x\ > 0, 
|Vuo| 2 =x 2 X2 on d{uo > 0} 

in the case x\ = x® — 0. 

Next, let us try to compute u : In the cases where > 0, the homogeneity of u 
and its harmonicity in {uq > 0} imply the following: if x 2 > 0, then each connected 
component of {uq > 0} is a half-plane passing trough the origin. If x 2 = 0, then 
the fact that uq must be harmonic in {^2 < 0}, implies that {uq > 0} is a cone 
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with vertex at the origin and of opening angle 120° symmetric with respect to and 
containing {(0, t) : t > 0}. 

In the cases where x\ = 0, solving the resulting ODE leads to hypergeometric 
functions and is slightly awkward, so we will instead use, in each section of the unit 
disk where uq > 0, the velocity potential (f> defined by 

d\4> = — d 2 u, d2(f> = d\u. 

X\ X\ 

In the case x\ > we obtain that <fi(p sin 9, pcosO) is homogeneous of degree 1 and 
is on the unit circle given by a linear combination of Pi(cos#) and $t(Q%(cos9)), 
where Pi and Qi are the Legendre functions. Now P\{x) = x and dlQi is a strictly 
convex function with singularities at —1 and 1, so that it is not possible that 

aP[{x) + (RQi)'(as) = aP[(y) + (HQi)'(y) for x ? y e (-1, 1). 

It follows that there can be at most one free surface point of the solution aP\ (cos 8)+ 
dlQi(cosd) in (0, 7r), but then the solution would have at least one singularity in 
the interval [0, tt]. Thus the only solution possible is a P\ (cos 9) — a cos 6, so that 
4>{x) = (JX2 and Uq(x) = cx\, where c and a are non-negative constants. The 
statement about the density follows as y = 1 m i u o > 0}. 

In the case x% = we obtain that ^(psin^jpcos^) is homogeneous of degree 
3/2 and is on the unit circle given by a liner combination of _P 3 / 2 (cos#) and 
P 3 / 2 (— cos9), where P 3 / 2 is the Legendre function. It is well known that P3/2 has 
only one singularity at —1 and that P3/2 has in (—1, 1) a unique zero z G (—1,0). 



By Lemma 3.9 we obtain as in the last case that aP 3 / 2 (cos9) + f3P 3 / 2 (— cos9) can 
have at most one free surface point in (0,7r). but then the solution would have at 
least one singularity in the interval [0, tt] unless /3 = 0. The fact that the singularity 
and the unique zero are both contained in [—1, 0) implies therefore that either 

cf>{p sin 6, p cos 6) — cr j o 3 / 2 P 3 / 2 (cos0) in {0 < 9 < arccos(z )} 



or 



4>(p sin 9, p cos 9) = <jp 3 / 2 P 3 / 2 (— cos 9) in {arccos(— z ) < 9 < tt}. 

However the free surface must not intersect {x 2 < 0}, so that we obtain that the 
only admissible solution is 

0(psin0, pcos9) = <jp z ^ 2 Pz/ 2 {— cos6>) in {arccos(— z ) < 9 < tt} 

for some nonzero constant a. Switching from the velocity potential back to uq we 
obtain the statement about uq as well as the density. 

Last, consider the situation when the set Bg(z) n {uq > 0} has two connected 
components. The computations of uq in the respective cases show that this is only 
possible for x\ > and x 2 > 0. The argument for (3.331 yields in this case that 
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the constant values of |Vmo| 2 on either side of 8{uq > 0} are equal. This concludes 
the proof. □ 



Lemma 3.11. Let u be a weak solution of (3.1 1 such that u = in {x 2 < 0} and 
suppose that 

IV^I 2 < + . 
< XixZ in \l. 

Xl 

Then x 2 — 0, x\ > and M X2 (0+) = imply that u = in some open 2- 
dimensional ball containing x° , while x\ = x 2 = M XlX2 (0+) = implies that 
u = in for some 5 > 0. 

Proof. Suppose towards a contradiction that a; e d{u > 0}, and let us take a 
blow-up sequence 

u m {x) := u(x° + r m x)/rK 2 

converging weakly in W^ 2 (R 2 ) to a blow-up limit uq in the case that x\ > 0, and 
a blow-up sequence 

u m (x) := u(x° + r m x)/r^l 2 

converging weakly in W^ 2 oc (T{? + ) to a blow-up limit u in the case that X® = 0. 
Proposition 3.10 shows that uq = in R 2 . Consequently, 

«- div (-g — Vu m )(S 2 ) > / y/x^dH 1 in the case a;? > 0, 

(3.37) 

<- div (— Vu m )(B%) > / y/x^dn 1 in the case x? = 

X l " iB+n3«d{llm>0} 

as m — y oo. (Recall that div (^V«) is a nonnegative Radon measure in f2.) On the 
other hand, there is at least one connected component V m of {it m > 0} touching 
the origin and containing by the maximum principle a point x m 6 dA, where 
A = (-1, 1) x (0, 1) in the case x? > and A = (0, 1) x (0, 1) in the case x? = 0. If 
max{x2 : x £ V m D dA} -/> as m — > oo, we immediately obtain a contradiction to 
(3.371. If max{x 2 : x G Kn n <9A} — > 0, we use the free-boundary condition as well 
as \Vu\ 2 /x 2 < x 2 to obtain 

= div(- 5 — Vu m )(V m nA) < [ y/x^dH 1 - f yfx~idU x 

X\ + T m Xi JV m ndA JAndrodVm 

in the case x^ > 0, 

= div (—Vu m )(V m n A) < / y/x^dn 1 - [ yfx^dU 1 
x i Jv m ndA JAnd Icd V m 

in the case x? = 0. 

However J v n9A y fx 2 d'H} is the unique minimiser of J gD ^Jx 2 dW 1 with respect to 
all open sets D with D — V m on dA. So V m cannot touch the origin, a contradiction. 

□ 
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Theorem 3.12 (Curve Case). Let u be a weak solution of (3.1) satisfying 

|Vu| 2 

< Ca^la^l locally in O, 

and let x° G Q be such that x^X® = 0. Suppose in addition that d{u > 0} fl B^~(x°) 
is in a neighborhood of x° a continuous injective curve a : I — > R 2 such that 
a = (<7i, (72) and tr(0) = X . Then the following hold: 
(ii) Stokes corner: If x® > 0, x% — and 



M X2 (0+)=X 1 / a; 2 X{(psine,pcose):-7r/3<e<ir/3}^) 



then (cf. Figure ^ o~i(t) ^ x\ in {— ti,ti)\{fS\ and, depending on the parametriz 
tion, either 

V ^ 1 A V g 2(t) 1 

lim — k = and lim — „ = =, 

t-yo+a^t) - x\ V3 t^o- a x {t) - x\ y/3' 



lim — , 2 S ^ — 7t = = and lim 



t-*)+ ai(t) -x\ ' t^o- ai (t) - x\ 




u = 



Figure 4. Stokes corner (x° > 0,^2 = 0) 
(i 2 ) Ifx\ > 0, x% = and M X2 (0+) = x\ f B + x% dx or M X2 (0+) = x\ J B + x^ dx, 
then (cf. Figure [5J) ^ x\ in (— 1±, t\) \ {0}, a\ — x\ changes sign at t = and 

lim „ = 0. 

o~ 1 (t) - arj 
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u > 




u = 



Figure 5. Horizontal point (x\ > 0, xE] = 0) 
(13) In £/ie case x\ > 0, x\ = and M X2 (0+) = — which is according to Lemma 



3.11 not possible at all provided that u = in {x 2 < 0} and the sharp Bernstein 



inequality holds — , then ai(i) 7^ x\ in (— ti,ti) \ {0}, o"i — x\ does not change its 
sign at t = 0, and 

ton g f) = 0. 

t->-0 CTl (t) - 

(mi) // ^2 > 0, x\ = and M X1 (0+) = x\ J g + x\ dx, then (cf. Figures 
02 (i) 7^ ^2 * n (0, ti) and 

ton g*W = o 

or cr 2 (i) 7^ in (— ii,ti) \ {0}, a 2 — x® changes sign at t = and 

lim °\® Q = 0. 

(M2) The case x\ > 0, x\ = and M a:i (0+) = is not possible. 
(mi) Garabedian corner: If x® = x\ = and 

M X12:2 (0+) = f xix 2 dx, 

J B+n{P^ /2 (-cos6)<0} 

then (cf. Figure [$[) o"i(i) 7^ in (0, ti) and, 



lira °" 2 ^| = tan(7r/2 — arccos(— zo))- 

i->0+ 0"i(t) 



(m'2) //a;? = x\ = and 



M Xi:C2 (0+)= / ni+diorM Xll2 (0+)= / ^a^da:, 
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Figure 6. Downward Vertical Cusp {x\ — 0,2;° > 0) 




FIGURE 7. Upward Vertical Cusp {x\ = 0, x\ > 0) 



then (cf. Figure 10) <J\{t) ^ in (0,ii) and 



lim ^44=0. 
t-yQ cj i (t) 
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FIGURE 8. Double Vertical Cusp (x° = 0, x° > 0) 





u = 




u > 



Figure 9. Garabedian corner {x\ 



0) 



(7n the subsequent sections of the present paper we will analyze the precise asymp- 
totics of the velocity field in the case M XlX2 (0+) = J B + X1X2 dx.) 



(m 3 j 



Ifx\ 



and M XlX2 (0+) = — which is according to Lemma 



3.11 not possible at all provided that u = in {x 2 < 0} and the sharp Bernstein 
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U > 




U = 



FIGURE 10. Horizontal point (x° = x° = 0) 

inequality holds — , then <Ji(t) ^ in (— t\,ti) \ {0}, and 

lim ^4=0. 

Remark 3.13. Although we omit a proof in the present paper, a perturbation of 
the frequency formula in |2T] (see [2D]) can be used to prove that, if x\ > and x® = 0, 
then case M X2 (0+) = x\ J B + x\ dx is not possible. Case (tit) seems possible as we 
have a nontrivial homogeneous solution. We do at present not have an existence 
proof for the cusps suggested here. 

Proof. We prove the claimed results only case (iii), when x\ = x\ = 0, the analysis 
in the other cases being similar. For each y = (2/1,2/2) € R 2 with y\ > and 
(j/i , 2/2 ) 7^ (0,0), we define 9(y) € [0,ir] by the relation 

G/1,2/2) = (p(y)sine(y),p(y) cos6(y)). 

We now consider the set 

£ = {#0 G [0, 7f] : there is t m — > such that 9(a{t m )) 9q as m —> 00}. 

Note that in fact C C [0,7r/2], since the free boundary d{u > 0} is contained in 
{x 2 > 0}. 

We now claim that: The set C is a subset o/{0, 9* , 7r/2} ; where 9* — arccos(— zq) 
is the angle corresponding to the Garabedian cone. 

Indeed, suppose towards a contradiction that a sequence ^ t m — » 0+,m — > 00 
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exists such that 9(a(t m j) — > 9 £ C \ {0, 9* , n/2}, let r m := \a(t m )\ and let 

u(r rn x) 



5/2 



For each p > such that B := ^(sin^o, cos#o) satisfies 

= Bn({(a,O) : a £ R+} U {(0,a) : a € R+} U {(asin#*,acos#*) : a € R+}) 



we infer from the formula for the unique blow-up limit u (see Theorem 3.101 that 
the convergence of measures 

(div (—Vu m ))(B) -> (div (— Vit ))(-B) = as m -> oo. 

On the other hand, 

div (—Vu m ) = ^Em 1 [d{u m > 0}, 

X\ 

which implies, since B n d{u m > 0} contains a curve of length at least 2p — o(l), 
that 

<— (div ( — \7u m ))(B) > c(8 ,p) as m —> oo, 

where c(6o,p) > 0, a contradiction. This proves the property claimed. 

Now, a continuity argument yields that £ is a connected set. Consequently the 
limit 

i = lim 9Mt)) 

exists and is contained in the set {0, 8* ,ir/2}. In what follows, we identify the value 
of £ in terms of the value of M XlX3 (0+). 

Suppose first that M XlX2 (0+) = J B + n ^ P , (- C os6)<o} XlX2 Then, by Propo- 
sition |3T0] the blow-up limit is 

Uo(psm9, pcos9) = r 3 Ut(6). 

Since (div (^-VMo))(-B 1 / 1O o(sin6'*,cos0*)) > 0, it follows that we cannot have £ £ 
{0, 7r/2}, and therefore we must have £ = 9* . This proves case (mi) of the Theorem. 

Suppose now that M XlX2 (0+) £ | f g + X\X^ dx, J B+ x\X^ Then the blow- 

up limit is uq(x) = 0. The same argument given earlier in the proof shows that 



£ 7^ 9*, so that necessarily £ £ {0,7r/2}. But then the formula in Lemma 3.8 that 
M XlX2 (0+) = lim r- 4 / Xl x 2 X{u>o}dx, 

shows that £ = implies M XlX2 (0+) = 0, while £ = it/2 implies that M XlX2 (0+) £ 
{ Ib + XlX t dx, J B + X1X2 dx, oj. However, the possibility that M XlX2 (0+) = and 
£ = is ruled out by the argument in the proof of Lemma [3.11| even in the absence 
of the strict Bernstein condition. This proves the cases (m'2) and (m'3) of the 
Theorem. 

□ 
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4. Frequency formula 



From now on we will focus on the case x\ = x 2 = 0, u = in {x 2 < 0} and 
M XlX2 (0+) = J B + x\x 2 dx, in which we will derive a precise asymptotic profile of 
the velocity. 



Theorem 4.1 (Frequency Formula). Let u be a variational solution of (3.1), and 
let 5 := dist(0, <3fi)/2. Let, for any r E (0, 5), 



D(r) 



Jb+(o) 7l\Vu\ 2 dx 



and 



V ^ = r 



IdB+(0) xi 



u 2 dU x 



Then the 'frequency" 

H(r) = D(r) - V(r) 

r Js+(o) (stI Vu I 2 + x i x 2(X{u>o} -1) 



dx 



/(5B+(0) xi 

satisfies for a.e. r € (0,6) the identities 
H'(r) 



L,,2 



u 2 dU l 



1 

r JdB + (0) x l 

2 



r(Vu • v) 

W(o) ir w2dH1 



1/2 



£>(r) 



1/2 



+ -V^r) + -V(r) ff(r) - - 



cm 1 
(4.1) 



and 
H'{r) 



1 



r J9i3 + (0) ^1 



r(Vu ■ i/) 



(/sb+(o) ir" 2 ^ 1 



1/2 



(W(o) k u2dnl 



1/2 



2 , , / , , 5 
+ -V(r) \B(r) - - 



Proof. Note that, for all r € (0,6), 



H(r) 



r 4 I(r) - j B+ xix 2 dx 



r- 5 J(r) 
Hence, for a.e. r £ (0,6), 

(r- A I(r))' (r^I(r) - f B + x x x 2 dx) (r- 5 J(r))' 



dH 1 
(4.2) 

(4.3) 



H'(r) = 



r 5 J(r) 



r 5 J(r) 



r 5 J(r) 
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Using the identities (3.191 and (3.121 with a — —5, we therefore obtain that, for 
a.e. r G (0,5), 



H'(r) = 



J+(0) u an 



. , , . v \ ( 2r hBt(o)k u{ ^ u - v ^ dnl ~ b hBt(Q)k u2dnl 
- (D(r) - V(r))-^ ^ a 

r Jc 

' L)(r) 



dB?(0) x 



-U 2 dU 1 



JdB+(0)xl u2dnl 



'(D(r)-V(r))(D(r) 



where we have also used the fact, which follows from (3.141, that 

r fdB+(o) £«(Vu ■ i/) dH 1 



D(r) = 



SdB+(o) k u2dV} 



(4.4) 



(4.5) 



Identity (4.1) now follows by merely rearranging (4.4 1, making use again of (4.5) 
and the fact that D(r) = V(r) + H(r). 



Since (4.1) holds, it follows by inspection that (4.2) holds if and only if 
1 



[r(Vu • v) - D(r)u] dU 1 + V 2 (r) 



1 

?+(0) x i 



vfdU 1 



{ — [r(\7u-v)-H(r)u] 2 d1-L l . 

JdB+(0) x l 



(4.6) 



However, (4.6 1 is easily verified as a consequence of (4.5) and the fact that D(r) 
H(r) + V(r). In conclusion, identity (4.2 1 also holds. 



□ 



Theorem 4.2. Let u be a variational solution of (3.1) such that u = in {X2 < 0}, 
let x° = 0, suppose that M XlX2 (0+) = J B +xixfdx, and let 5 := dist(0, dQ)/2. 
Then the following hold: 

(i) H(r)>\ for all r G (0,5). 

(ii) The function r i— )> r~°J(r) is nondecreasing on (0,5). 

(Hi) The function H is nondecreasing on (0,5), and has a right limit H(0+), 
where H(0+) > 5/2. 

(iv) r ^ W 2 (r) G V-%5). 



Proof, (i) The monotonicity, which follows from Theorem 3.2 of the function M XlX2 
ensures that, for all r G (0, 5), 

5 



r -*j( r ) _ ~-r~ b J(r) > / Xl x+ dx. 



(4.7) 



Using (4.3 1, the above inequality may be rearranged in the form of the claimed 
result. 
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(ii) Plugging a = — 5 into (3.121, using also (3.14), and then (4.7), we obtain, 
for a.e. r £ (0, S), 



(r" 5 J(r))' = - I r" 



> 2r~ 



|Vu| dx 



1 



B+(Q) x i' ' 2 JdB+(Q) x l 

x\x 2 {l ~ X{u>o})dx > 



u 2 dU x 



B+(0) 

which implies the claimed result. 

(iii) The monotonicity of H on (0,8) is a consequence of (4.1l and (i). The 
remaining part of the claim is immediate. 



(iv) The claimed result follows from (4.1l and (iii) 



□ 



5. Blow-up limits 
The Frequency Formula allows passing to blow-up limits. 

Proposition 5.1. Let u be a variational solution of ( |3.1| such that u = ir, 
{x2 < 0}, let x° = 0, and suppose that M XlX2 (0+) = J B + x\x\ dx. Then: 

(i) There exist lim r _ ! .o+ V(r) = and lim r _ ) .o+ D(r) = H(0+). 

(ii) For any sequence r m —> 0+ as m —} oo, the sequence 

u(r m x) 



v m (x) :-- 



(5.1) 



bounded in W^ 2 {BJ) 



(iii) For any sequence r m — > 0+ as m — > oo such that the sequence v m in (5.1) 
converges weakly in W^ 2 (B^) to a blow-up limit vq, the function vq is homogeneous 
of degree H(0+) in , and satisfies 

v > in B lt v Q = in Bf n {x 2 < 0} and [ —v 2 dU x = 1. 

JdB+ Xl 

Proof. We first prove that, for any sequence r m — > 0+, the sequence v m defined in 



(5.1 1 satisfies, for every < r < a < 1 
1 



B+\B+ Xi 



\x\ [Vv m (x) ■ x - H(0+)v m (x)] dx -t asm-4oo. (5.2) 



Indeed, for any such r and cr, it follows by scaling from (4.2) that, for every m such 
that r m < 5, 



2 [ 1 

r JdB+ x i 



r(Vu m • v) 



— ~H(r m r)- 



< H(r m o~) — H(r rn T) — > as m — > oo, 



1/2 



<m}dr 
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as a consequence of Theorem 4.2 (iii). The above implies that 



1 

j+ x i 



r(Vv ri 



as to — > oo 



1/2 



H(0- 



(X 



1/2 



<m}dr 



(5.3) 



Now note that, for every re (r, <r) C (0,1) and all m as before, it follows by using 



Theorem 4.2 (ii), that 



3+ %i 



-v 2 dH 1 



>0B+ x x 



^U 2 dH 1 



Id 



< r 5 < 1. 



Therefore (5.2 1 follows from (5.3), which proves our claim. Let us also recall (4.5 1 



We can now prove all parts of the Proposition. 

(i) Suppose towards a contradiction that (i) is not true. Let s m — > be such 
that the sequence V(s m ) is bounded away from 0. It is a consequence of Theorem 
OJiv) that 

V{r) — > as m — > oo. 



mm 

r£[s m ,2s„ 



Let t m £ [s m , 2s m ] be such that V(t m ) — > as m — > oo. For the choice r m :— t m for 
every m, the sequence v m given by (5.1) satisfies (5.2 1. The fact that V(r m ) —> 



implies that D(r m ) is bounded, and hence that v m is bounded in W^ 2 (B^). Let vq 
be any weak limit of v m along a subsequence. Note that by the compact embedding 
W^' 2 (B^~) ^-s- L 2 (dB^), vq has norm 1 on L^dB^), since this is true for v m for 



all 777. It follows from (5.2 1 that vq is homogeneous of degree H(0+). Note that, 
by using Theorem |4.2| (ii), 

%* 4 Js+ ^1^2(1 - X{»o}) dx 



V(s m ) 



±u 2 dH 1 

Xl 



< 



n 4 J B + X i x 2^ - X{u>a}) dx 



< 



1 Job? 



»/2 

±u 2 dU 1 

Xl 



2 X 



f)B^ 



^rU 2 dU 



9B+ /2 



Since, at least along a subsequence, 

—vLdH 1 

9B+2 Xl 



T V(r m ). 



(5.4) 



j+ Xl 

'1/2 



-v 2 dH 1 > 0, 



(5.4 1 leads to a contradiction. It follows that indeed V(r) — ► as r — >• 0+. This 
implies that D(r) -> #(0+). 
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(ii) Let r m be an arbitrary sequence with r m — > 0+. The boundedness of the 
sequence v m in W^' (Bi) is equivalent to the boundedness of D(r m ), which is true 

by (i). 

(iii) Let r m —> 0+ be an arbitrary sequence such that v m converges weakly to 
vq. The homogeneity degree H(0+) of vq follows directly from (5.2). The fact 
that J dB + -^Vq dH 1 = 1 is a consequence of J dB + ^-"J, dH 1 = 1 for all m, and the 
remaining claims of the Proposition are obvious. The homogeneity of vo, together 
with the fact that vo belongs to W^ 2 (B^), imply (in two dimensions) that vo is 
continuous. 

□ 



6. CONCENTRATION COMPACTNESS 



In the present section we will prove a concentration compactness result which 
allows us to preserve variational solutions in the blow-up limit at degenerate points 
and excludes concentration. In order to do so we combine the concentration com- 
pactness result of J.-M. Delort |S] with information gained by our Frequency For- 
mula. In addition, we obtain strong convergence of our blow-up sequence which is 
necessary in order to prove our main theorems. 



Theorem 6.1. Let u be a variational solution o/(3.1l such that u = in X2 < 
and M XlX2 (0+) = J B+ x x x^ dx. Let r n 



— > 0+ be such that the sequence v m given 
by (5.1) converges weakly to vo in W^' 2 (B^~). Then v m converges to vq strongly in 
^i/loc(-®lA{0})> vq is continuous on and div (^Vuq) * s a nonnegative Radon 
measure satisfying uodiv (^V«o) = in the sense of Radon measures in B^~ . 



Proof. Note first that the homogeneity of vo given by Proposition |5.1[ together with 
the fact that vo belongs to W^ 2 (B^), imply that vq is continuous. 

Let a and r with 0<T<er<lbe arbitrary. We know that div (-^-Vv m ) > 
and div (^Vw m )(Bj. +1 y 2 ) < C\ for all m. We regularize each v m to 



wher 



and 



is a standard mollifier such that 

1 f 1 

div ( — Vv m ) > 0, / div ( — Vv m ) <C 2 < +oo for 
xi J B + xi 



all m, 



1 



(Vv m - Vv m ) 



LHB+) 



as to — > oo. 



Let us now consider the velocity field in three dimensions 
V m (X,Y,Z) := 



— cbwmCosi?, d2V m sm , d, — div n 

X\ X\ X\ 
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where (X, Y, Z) — (x\ cos x\ sin $, X2), as well as their weak limit 
V(X,Y,Z) := ( - — d 2 icostf,- — d 2 vsmtf, — div 

\ X\ X\ X\ 

We have that V m is divergence free and satisfies 

curl V m = co m = (-sintf,costf,0)a m in B 2 (0) 
with a non-negative function a m that is bounded in L l {B a ). It follows that 

V? = A m \d z uj? 

where A~ \ is the inverse of the three dimensional Laplace operator with averaged 
Dirichlet boundary data V" 1 , more precisely 

l + cr 1 + CT 

A-i/= T ^/ 2 f G R fdxdR+-^— f 2 f V?VG R -vdH\ 

1 - J a J Br I - (J J a J 9Br 

where Gr is Green's function with respect to the Laplace operator in Br. From 
the proof of Proposition 3.2], where [HI (3.6)] holds with v\ replaced by V™ and 
wf replaced by w™ but the remainder terms wf given by Greens formula in B a . we 
infer that 

VrV 3 m - VM weakly in I%jB a ), 
V 2 m V 3 m - V 2 V 3 weakly in L 2 oc (B a ), 
(Vrf + (V 2 m f ~ (V 3 m f ~* (Vi) 2 + {V 2 f - (V 3 ) 2 weakly in L 2 oc (B a ); 

note that as in |8] the remainder terms converge strongly in L 2 oc {B a ). 
It follows that 

—div m d 2 v m -> —div Q d 2 v (6.1) 

X\ X\ 



and 



- ((d!v m ) 2 - (d 2 v m ) 2 ) -> — ((d lVo ) 2 - (d 2 v ) 2 ) 

X\ ' X\ 



in the sense of distributions on B^ as m — > 00. Let us remark that in contrast 
to the true two-dimensional problem, this alone would not allow us to pass to the 
limit in the domain variation formula for v m ! 



Observe now that (5.2 1 shows that 

Vv m (x) ■ x - H(0+)v m (x) -> 
strongly in L 2 W {B+ \ £?+) as m — > 00. It follows that 

dxVjnXx + d 2 v m x 2 -> diVQXx + d 2 v x 2 
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strongly in L 2 W {B+ \ B+) as m — > oo. But then 



1 

B+\B+ x l 



+ div m d 2 v m x 2 )ri dx 



1 



{diVodiv a xi + div d 2 v X2)r] dx 



for each rj £ Cq(B+ \ B T ) asm-) oo. Using (6.1 I, we obtain that 



:\B+ 



{div m ) 2 r]dx 



{divof'qdx 



for each < rj £ Cq(B+ \ B T ) as to — > oo. Using once more (6.1 1 yields that Vw m 
converges strongly in l? w loc (B+ \ B^). Since a and r with < r < <t < 1 were 
arbitrary, it follows that Vi> m converges to V«o strongly in L 2 W ^(-B^ \ {0}). 
As a consequence of the strong convergence, we see that 

/ — V(r?u Q ) • Vv = for all 77 e ^(5+ \ {0}). 

Combined with the fact that vq = in Bf n {0:2 < 0}, this proves that «oAwo = 
in the sense of Radon measures on B^~ . 

□ 



7. Degenerate points 



Theorem 7.1. Let u be a weak solution of ( |3.1[ ) suc/i £/ia£ u = in X2 < and 

M XlX2 (Q+) = J B + X1X2 dx, let the free boundary d{u > 0} n B^ be a continuous 
injective curve a — io\,a^) such that er(0) = 0. Then ai(t) ^ in [0,ti) \ {0}, 



lim 
t->o a x (t) 



and 



u(rx) 



x\x 2 



95+ (0) 



u 2 dn l 



\J fdB + (0) x \ x \ dH 1 



strongly in W^ ioc {B x \ {0}) and weakly in W ■ (B x ). Moreover 



u(rx) 
u(rx) 



-> in L 2 W (B+) for a G (0, 2) and 
is unbounded in L 2 U (B^) for a > 2. 



Proof. Let r m — > 0+ be an arbitrary sequence such that the sequence v m given 
by (5.1 1 converges weakly in W^ 2 {B X ) to a limit vq. By Proposition 5.1 (iii) and 



Theorem 6.1 vq ^ 0, Vq is homogeneous of degree H(0+) > 5/2, vq is continuous, 
vo > and vo = on {xi = 0} and in {x 2 < 0}, wodiv (;^-Vt>o) = in as 
a Radon measure, and the convergence of v m to vq is strong in W^ 2 oc {B x \ {0}). 
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Moreover, the strong convergence of v m and the fact proved in Proposition 5.1 (i) 
that V(r m ) — V as tyl — v oo imply that 



0= / (— |Vu | 2 div <j>- 2\7u D(j>Vvo 

J R 2 \Xl 



for every cf> g Cq({xi > 0} fl {X2 > 0}; R 2 ), so that even an analysis in the case of 
{u = 0} consisting of infinitely many disconnected components (similar to that in 
|21| ) would be possible in principle. However the structure here is more complicated. 
For that reason we confine ourselves to the assumed injective curve case. 



As in the proof of Proposition 3.10 we will use in each section of the unit disk 



where vq > the velocity potential (f> defined by 

d X (j) = —d 2 v , d 2 (f> = — — div . 

Xi Xi 

We obtain that 4>(p sin 6, pcosB) is homogeneous of degree m = H(Q+) > 5/2 
and is on the unit circle given by a linear combination /(cos#) = aP m (cos9) + 
l3P m {— cos#), in the case that the Legendre function P m and the function P m (—x) 
are linearly independent, and /(cos 0) — aP m (cos 0) +(3di(Q m (cos 9)) in the case the 
Legendre function P m and the function P m {—x) are linearly dependent. Moreover 
(1, 0) is a free boundary point of vo so that /'(0) = 0, which implies a — f3 in the 
case of linear independence. 

(ii) implies that for any ball B CC n 



On the other hand, Theorem 

u{rx) 



3.12 



{x2 > 0}, v r = , u ( rx ) > in B. Consequently div ( — Vwo) = in 

V /r l JeB+Co) u2<iH1 

{x\ > 0} fl {x2 > 0}. However, if there is a free boundary point x in (0, 1) x (0, 1) 
then by homogeneity the half line connecting that point to the origin consists of free 
boundary points, so that (div Vvo)){Bs(x)) > for each 5 > 0, a contradiction. 
Thus ctP' m + PQ' m must be either strictly positive or strictly negative in (0, 1). 

In the case /(cos#) = a(P m (cos 8) +P m (— cos0)) we obtain now a contradiction 
to the fact that P m is bounded at 1 and has a singularity at —1. 

In the case that P m is an even function, we obtain from P' m (Qi) — mP m -i(0) = 

m V™ ari A r>' <f\\ — mO , — »»r 3/2 tan(7r(m-l)/2) 

T{^)T(^+V, and - mt-im-lKU) - -( m _ 1)r( 2^Hi )r (ni^i) 

(see| http://functions.wolfram.com/07.07.20.0006.01[ 
http: 



http 
|http 



//functions.wolfram.com/07.07.03. 0001.01, 
functions.wolfram.com/07.10. 20. 0003.01, 



//functions.wolfram.com/07.10. 03. 0001.01) 



that m is an even integer > 2 and that /3 = so that / is up to a nonzero 
multiplicative constant the Legendre polynomial P m . But, using [3J Corollary on 
p. 114] there is only one even integer > 2 such that P m has no critical point in 
(0, 1), namely m = 2. We obtain f(x) — c 2 P 2 (^) = c 2 |(3a; 2 - 1). 
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.",5 



In order to obtain the claimed growth we calculate for u r {x) = u(rx)/r a and 



a.e. r G (0,S), using (3.14 1, 
1 



-ul dU 1 = 



9B+(0) x l 



— |Vu r | 2 dx — a 



> * f. 

— r J t 



1 „,2 



uidW 



dB+(0) xi "r 
r JdB+(0) xi w r 



B+(0) x l 

a 6(0,2), 



1 

dB+(0) x l 



uldU 1 



Integrating we obtain the result. 



□ 
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